How to observe dipolar effects in spinor Bose-Einstein condensates 
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We propose an experiment which proves the possibility of spinning gaseous media via dipolar 
interactions in the spirit of the famous Einstein-de Haas effect for ferromagnets. The main idea is 
to utilize resonances we find in spinor condensates of alkali atoms while these systems are placed 
in an oscillating magnetic field. A significant transfer of angular momentum from spin to motional 
degrees of freedom observed on resonance is a spectacular manifestation of dipolar effects in spinor 
condensates. 



Historically, first attempts at proving the relationship 
between magnetism and angular momentum (i.e., the 
Ampere's hypothesis on "molecular currents") were al- 
ready performed in the 19th century [l[. They were all 
unsuccessful because of difficulties in measuring the tiny 
forces involved in the process. The efforts continued in 
the beginning of 20th century resulted in the emergence 
of a new class of effects named as magnetomechanical. 
The most known is certainly the Einstein-de Haas effect 
0] in which a magnetized ferromagnetic rod is forced to 
rotate when its magnetization is reversed. In the exper- 
iment the magnetomechanical ratio is measured and its 
value can be directly related to the motion of electrons in 
the rod. It tells us how large a portion of magnetization 
comes from the spin of electrons and how large from the 
orbital motion of electrons. 

Experiments with neutral "molecular currents" would 
be of particular interest in proving the direct and sole 
relation between the spin and the magnetization of the 
system. The recent realization of chromium condensates 
[3|] has launched huge interest in ultracold dipolar sys- 
tems 0, [B| which would be an ideal candidate for such 
studies. Although spectacular features due to dipolar 
forces related to expansion Q and collapse @] were al- 
ready observed, the Einstein-de Haas effect so far remains 
elusive. A route towards observing the Einstein-de Haas 
effect in chromium condensates via controlling the mag- 
netic field was recently suggested [H, Q and first experi- 
ments demonstrated the possibility to control the dipolar 
relaxation rate in this way [Icj . 

It is tempting to test these concepts also in alkali sys- 
tems which represent the majority of experiments. It is 
certainly less obvious that alkali atoms (whose magnetic 
dipole moment is an order of magnitude lower than that 
of chromium atoms) might be a good candidate to ob- 
serve how spin is transmitted into orbital motion. How- 
ever, some authors suggested that the magnetic dipolar 
interactions could already lead to observable effects in 
condensates of alkali atoms 0, [lll - [l7| . Also a first exper- 
iment showed that a spin-1 87 Rb spinor condensate can 
exhibit dipolar properties [l8j |. 

In this Letter we demonstrate that controlling the 
dipolar interactions by using an oscillating magnetic field 



is a perfect way to observe the Eintein-de Haas effect in 
alkali condensates provided it is done under a resonance 
condition. In the following we explain our reasoning in 
detail. 

The equation of motion for a spinor condensate in the 
F = 1 hyperfine state in the mean-field approximation 
reads 

ih^(r) = ['H sp +H c + U d )i>(r), (1) 

where ^(r) = (^i(r), i/jq(t), V>-i(r)) T is a condensate 
spinor wavefunction and the effective Hamiltonian is split 
into three terms according to their origin. And so, the 
single-particle contribution is given by 



H sp — H sp -fiBF 



(2) 



where H sv = -|-V 2 

^ 2m 



V tr , P = (F x ,F y ,F z ) are stan- 
dard F = 1 spin matrices, B is an external magnetic field 
(assumed as directed along the z-axis), and [i is the mag- 
netic moment of an atom. The short-range interactions 
between atoms lead to the second term 

H c = c ^(r) V(r) + ^(^(r) F ^(r)) F , (3) 

where cq = 4irh 2 (2a,2 + ao)/3m and C2 = 4-7r?i 2 (a2 — 
ao)/3m determine the strength of spin-preserving and 
spin-changing collisions, respectively [l9| : do (02) is the 
scattering length of colliding atoms for the channel of 
total spin equal to zero (two). Finally, the contribution 
corresponding to the long-range dipole-dipole interaction 
takes the form 



d 3 r'ip^(r')V d (r-r')^(r'). 



(4) 



Here, the dipolar interactions V d can be written in the 
way of an expansion in the spherical harmonics functions 
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where Y2\(r) (with f denoting a unit vector in the di- 
rection of relative position of two interacting atoms) is a 
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spherical harmonics of rank-2 and £2. a defined as 



J 2,0 



-(^(r)F 2 (r')-F(r).F(r')/3) 



S 3l ±i - ±-(F z (v)F ± (v') + F ± (v)F z (v')) 



J 2,±2 



-F±(r)F ± (r>) 



(6) 



is rank-2 spherical tensor built of atomic spin operators 
(here, F±=F x ±iF y ). 

It follows from ([6]) that when two atoms collide the to- 
tal spin projection can change at most by 2 whereas the 
spin projection of individual atoms changes maximally 
by 1. Therefore, the atom can not be transferred directly 
from the nip = I to mp = — 1 component. The last row 
in ([5]) shows the processes in which the spin of each atom 
changes by +1 or —1. In addition, there are atomic col- 
lisions in which only one atom is transferred to the other 
Zeeman state (see the middle row in ([6|)) or the collisions 
that do not change the total spin projection (the first 
row in ©)• In the latter case, however, still the process 
transferring one atom from mp = to mp = 1 and the 
other one from mp ~ Q to mp = —1 states is allowed. 
According to ((SJ) in each case the change of the total spin 
projection is accompanied by an appropriate change of a 
relative orbital angular momentum of colliding atoms. 

In order to gain some analytic insights into the physics, 
we first assume that the transfer to the mp = — 1 state is 
small such that the component can be neglected. 
This assumption is justified by the full scale numeri- 
cal calculations below and allows to write the evolu- 
tion of the system as to be governed effectively by the 
2x2 Hamiltonian matrix. Further, by making a trans- 
formation into the interaction picture with the help of 
exp (-| J * H (t') dt'), where H = H' sp - \ p, B(r, t) with 
HL 



H sp + (cq + c 2 ) \ipi\ 2 + c \ip \ 2 and 



B(r,t) = B(t) + 6B(r,t) 
SB(r,t) = -(c 2 \^ Q \ 2 +n d n)/p 
one eventually obtains the Hamiltonian in the form 



(7) 



(8) 



provided the matrix elements Hiq = Hdw an d Hu = 
i uB(r,t) change slowly in space. 

The Hamiltonian <JSj) describes a two-level system as 
considered, for instance, in Ref. [2lj. Both Hiq and Hu 
are position-dependent functions however, which means 
that the spinor condensate evolving according to Eq. (TTJ) 
has been reduced, in fact, to the set of coupled two-level 
systems. The coupling comes through the off-diagonal el- 
ements, Hio, as well as through the diagonal ones, Hu, 
via Tidu- Looking for resonances we consider the oscil- 
lating magnetic field in the form (see Ref. (2l| ) 



Therefore, B(r,t) — Bo(r,t) + A cos uit, where Bo{v,t) — 
Bq+5B{y, t) can be interpreted as the position- and time- 
dependent nonoscillating part of the magnetic field. The 
correction SB(r,t) to the magnetic field originates from 
the contact interactions (the first term in 0) as well 
as the dipolar ones (the second term) which is the z- 
component of the magnetic field generated by the dipole 
moments of all the atoms outside the cloud of atoms. 

In the simplest (but not realistic) case H\o = const 
and SB(r,t) = and the spinor condensate effectively 
reduces to a single two-level system. Then going to the 
interaction picture with the coupling term treated as a 
perturbation (see Ref. [2l| for details) one finds that 
the evolution is governed by the Hamiltonian whose only 
nonzero elements are off-diagonal ones 
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(10) 



where J n (x) are the Bessel functions of the first kind. 
In a rotating wave approximation (RWA) all rapidly os- 
cillating terms in the sum in Eq. (fTU|) are neglected in 
comparison with the one, least oscillating, term. The 
resonance condition is then identified by 



nhw + fiB = 



(11) 



B(t) = B + Acoswt . 



(9) 



Within RWA and on resonance the population of mp = 
state oscillates with the frequency 257 such that Ml = 
\H\o\ I J n (uA/hu>)\. Therefore, to observe a significant 
transfer to the other state it is required to keep the fre- 
quency f2 as high as possible. This can be achieved by 
working in a strong-driving regime A ~ Bo (as opposed 
to weak-driving Rabi oscillations limit, A Bo) because 
then the Bessel function J n (uA/huj) ~ J n (pBo/huj) = 
J n (n) gets maximal values. 

To learn the values of parameters characterizing the 
oscillating magnetic field we use the formula (fTTj) for the 
static field Bo = lOmG which is reasonably large from 
experimental point of view. For a rubidium and sodium 
condensates one has a = /is/2. Assuming uj — 2n x 1Hz 
one gets n as 7000. To decrease the order of the resonance 
to n ~ 1 while keeping the same value of the static field 
one would have to increase the frequency u> by three or- 
ders of magnitude. Later, we will give arguments against 
increasing the frequency of an oscillating magnetic field. 

The spinor condensate is, however, equivalent to a set 
of inherently coupled two-level systems. Initially, all two- 
level systems are on resonance (condition (|11| is fulfilled) . 
Therefore, due to dipolar interactions, transfer of atoms 
to the mp = state begins. While this transfer occurs, 
the two-level systems get immediately out of resonance. 
They go off resonance in different ways, i.e., the nonoscil- 
lating part of the field for each two-level system depar- 
tures from Bo by different amounts. Fortunately, there 
exist side resonances according to (n ± k)hio + uB' = 



3 



where n and u> satisfy the condition ([TTjl for a given mag- 
netic field Bq and B' = Bo^fkhco/ fi (k is an integer num- 
ber). Therefore, further transfer to the mp = state can 
be sustained assuming the resonances are broad enough 
so they mutually overlap. The single two-level model just 
discussed gives answer to what influences the width of the 
resonances. 

Let's departure from the resonance condition Eq. (fTTj) 
by changing the static magnetic field: Bq — > Bq + 8 Bo. 
Then the least oscillating term in the off-diagonal element 
is given by 



Hiq J n 



/.iA 



-i^i 5Bq t/h 



(12) 



and the solution of an equation of motion for a two-level 
system with (fT2|) as the off-diagonal term reads 



if2_ t 



-iQ + t\ 



(13) 



where C = Hiq J n (fiA/huj)/h, w = (iSBq/H, f2 = 
■\/|C| 2 + (w/2) 2 , and fl± = f2±£j/2. Hence, the maximal 
off resonance transfer is calculated as 



\Mt) 



2 

max 



1 



1 



2|Hio| J„(M/fiu) 



2 ■ 



(14) 



First, it is clear from (|14p that the maximal transfer 
becomes larger when the off-diagonal element |iiio| in- 
creases, i.e., when the maximal atomic density increases. 
Second, it occurs from (fT4| that the transfer gets higher 
when the frequency of an oscillating magnetic field is low- 
ered while the static field Bq is kept constant (since the 
Bessel functions behave as nJ n (n) — > oo when n — > oo). 
Therefore, in both these cases the resonances broaden 
and for large enough \Hiq\ and/or small enough uj the 
resonances start to overlap significantly supporting in 
this way the transfer of atoms to mp = component. 
However, one has to remember that the formula (| 14[) is 
derived assuming the presence of only one term in the ex- 
pansion (fT0|) . Considering the full equation of motion for 
a two-level system (all terms included in the expansion 
(fTTJ)) ) it can be shown that, in fact, the transfer stops for 
small enough frequencies [22j . 

In the following we will demonstrate that the discus- 
sion of a two-level system presented above is indeed help- 
ful while looking for dipolar resonances by solving nu- 
merically Eq. |T]) (which is, generally, a search in three- 
dimensional parameter space). We prepare initially a 
system of N = 2 x 10 5 87 Rb atoms in a spherically sym- 
metric harmonic trap with frequency u>t r — 27r x 100Hz 
and in mp = 1 Zeeman state. We turn on the external 
time-dependent magnetic field with the static part and 
the amplitude of an oscillating term being of the order of 
lmG. The main frame in Fig. [T]shows the case of n = 645 
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FIG. 1. (color online). Relative transfer of Rb atoms to the 
t»F = and rriF = — 1 states as a function of time. Initial 
number of atoms in the tjif — +1 component is 2 x 10 5 . In 
the spherically symmetric trap with uj t r = 2n x 100Hz it cor- 
responds to the maximal density equal to 3 x 10 14 cm~ 3 . Pa- 
rameters of an oscillating magnetic field are: -Bo = — 0.92mG, 
A — 0.93mG, and uj/uitr = 1/100. Inset shows transfer of 
atoms for other resonance: Bo = — 1.43mG, A — 1.85mG, 
and lo /totr = 1/30. 



resonance whereas in the inset the order of the resonance 
equals n = 300. For both resonances the population of 
mp = — 1 component is negligible for shown time du- 
ration, hence the two-level system approximation works 
well in these cases. Note that the population of = 
state exhibits characteristic step-like behavior (see Ref. 
(2l| ) with jumps in the population occurring with the fre- 
quency of an oscillating magnetic field. The transfer to 
mp = state is caused purely by dipolar forces since the 
projection of orbital angular momentum per atom in this 
state equals h (in fact, a singly charged vortex is created 
in mp = component). 

Since the three-body loss rate constant for 23 Na atoms 
is a factor of 5 smaller than for 87 Rb atoms (2^1, the 
sodium spinor condensate seems to be more appropriate 
for the observation of dipolar resonances than the rubid- 
ium one. Lower losses allow for the higher atomic densi- 
ties (in fact, up to the order of 10 15 cm -3 ) Therefore, 
according to the two-level model higher transfers are ex- 
pected. In other words, a significant transfer for larger 
magnetic fields is still expected. In Fig. [5] we show the 
resonance for a sodium condensate at magnetic fields of 
about 3.5 mG. Indeed, the transfer is big and, what is 
also important, the width of the resonance is of the order 
of milligaus. Therefore, the quality of this resonance is 
much higher than those discussed in Refs. (il [l^. 

The quality of the resonance can be also improved by 
applying an inhomogeneous magnetic field as a nonoscil- 
lating part of external field. Surprisingly and opposite to 
what is predicted for the amplification of matter waves 
process in dipolar spinor Bose-Einstein condensates (see 
Ref. [25|]) the nonzero field gradient might enhance the 
population oi mp = state. It happens because the 
correction SB (see ((7|)) to the nonoscillating part of the 
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FIG. 2. (color online). Maximal (within 2.5s) transfer of 
23 Na atoms to the ttlf — state as a function of the static 
magnetic field |Bo|- Initial number of atoms in the tuf — +1 
component is 10 8 . In the spherically symmetric trap with 
u)tr = 2ir x 100Hz it corresponds to the maximal density 
equal to 1.08 x 10 15 cm -3 . Frequency and amplitude of an 
oscillating magnetic field are ui/uitr = 1/100 and A = 3.5mG, 
respectively. Inset shows transfer for Bq — — 3.4mG. 
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FIG. 3. (color online). Nonoscillating part of the magnetic 
field (as in 0) along the line shifted by 10.5 fim with respect 
to z-axis. The static magnetic field at the center of the trap 
equals Bo = — 14.0mG and its gradient is 2.8mG/cm. Other 
parameters are: ui/ujtr = 1/3 and A = 14.3mG. Successive 
curves (from top to bottom) correspond to times 0.2s, 0.4s, 
and 0.8s as marked by bullets in the inset. 

magnetic field, which is negative for sodium atoms, can- 
cels the increase of the field for some positive z. Even 
more, since the field gradient is present, still new regions 
in space serve as a source of atoms to go from mp = 1 
to mp — state (the n = 295 resonance is moving out of 
the center, see Fig. [3]). Therefore, the initial resonance 
is still active in regions in space somewhere at positive 
z. This resonance boosts transfer of atoms to vtif = 
state and at some time ignites the resonance of the order 
of n + 1 and even of n + 2 leading to further transfer 
(see inset in Fig. [3]). In fact, inhomogeneous magnetic 
fields support the transition of atoms in the case when 
this transfer is already stopped for uniform fields. 

In conclusion, we have shown how to control weak 
dipolar interactions with the help of an oscillating mag- 



netic field. Working in the resonant regime, which is 
experimentally accessible, one is able to transfer a sig- 
nificant number of atoms from the initial nonrotating 
state to the state with nonzero angular orbital momen- 
tum. Hence, the way to observing the Einstein-de Haas 
effect even in systems with very weak dipolar forces like 
ultracold alkali gases seems to be open. In fact, utiliz- 
ing the resonances opens new possibilities in studying a 
broad range of dipolar physics in alkali condensates. 
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